Intended for Phys. Rev. D 

Authors: Y. P. Guo, C. Z. Yuan, C. X. Yu 

Committee: M. Maggiora (Chair), F. Liu, X. T. Huan; 



Search for hadronic transition XcJ — > f]c^ f an d observation of 

XcJ ~+ KK7T7T7T 

M. Ablikim 1 , M. N. Achasov 5 , D. J. Ambrose 39 , F. F. An 1 , Q. An 40 , Z. H. An 1 , J. Z. Bai 1 , 
Y. Ban 27 , J. Becker 2 , M. Bertani 18A , J. M. Bian 38 , E. Boger 20 ' a , 0. Bondarenko 21 , I. Boyko 20 , 
R. A. Briere 3 , V. Bytev 20 , X. Cai 1 , 0. Cakir 35A , A. Calcaterra 18A , G. F. Cao 1 , S. A. Cetin 355 , 
J. F. Chang 1 , G. Chelkov 20 ' , G. Chen 1 , H. S. Chen 1 , J. C. Chen 1 , M. L. Chen 1 , S. J. Chen 25 , 

Y. B. Chen 1 , H. P. Cheng 14 , Y. P. Chu 1 , D. Cronin-Hennessy 38 , H. L. Dai 1 , J. P. Dai 1 , 
D. Dedovich 20 , Z. Y. Deng 1 , A. Denig 19 , I. Denysenko 20 ' 6 , M. Destefanis 43A ' 43C , W. M. Ding 29 , 
Y. Ding 23 , L. Y. Dong 1 , M. Y. Dong 1 , S. X. Du 46 , J. Fang 1 , S. S. Fang 1 , L. Fava 43S ' 43C , 
F. Feldbauer 2 , C. Q. Feng 40 , R. B. Ferroli 18A , C. D. Fu\ J. L. Fu 25 , Y. Gao 34 , C. Geng 40 , 
K. Goetzen 7 , W. X. Gong 1 , W. Gradl 19 , M. Greco 43A ' 43C , M. H. Gu 1 , Y. T. Gu 9 , Y. H. Guan 6 , 

A. Q. Guo 26 , L. B. Guo 24 , Y.P Guo 26 , Y. L. Han 1 , F. A. Harris 37 , K. L. He 1 , M. He 1 , 

Z. Y. He 26 , T. Held 2 , Y. K. Heng 1 , Z. L. Hou 1 , H. M. Hu 1 , J. F. Hu 6 , T. Hu 1 , G. M. Huang 15 , 
J. S. Huang 12 , X. T. Huang 29 , Y. P. Huang 1 , T. Hussain 42 , C. S. Ji 40 , Q. Ji 1 , X. B. Ji 1 , 
X. L. Ji 1 , L. L. Jiang 1 , X. S. Jiang 1 , J. B. Jiao 29 , Z. Jiao 14 , D. P. Jin 1 , S. Jin 1 , F. F. Jing 34 , 
N. Kalantar-Nayestanaki 21 , M. Kavatsyuk 21 , W. Kuehn 36 , W. Lai 1 , J. S. Lange 36 , C. H. Li 1 , 
Cheng Li 40 , Cui Li 40 , D. M. Li 46 , F. Li 1 , G. Li 1 , H. B. Li 1 , J. C. Li 1 , K. Li 10 , Lei Li 1 , Q. J. Li 1 , 
S. L. Li 1 , W. D. Li 1 , W. G. Li 1 , X. L. Li 29 , X. N. Li 1 , X. Q. Li 26 , X. R. Li 28 , Z. B. Li 33 , 
H. Liang 40 , Y. F. Liang 31 , Y. T. Liang 36 , G. R. Liao 34 , X. T. Liao 1 , B. J. Liu 1 , C. L. Liu 3 , 
C. X. Liu 1 , C. Y. Liu 1 , F. H. Liu 30 , Fang Liu 1 , Feng Liu 15 , H. Liu 1 , H. B. Liu 6 , H. H. Liu 13 , 
H. M. Liu 1 , H. W. Liu 1 , J. P. Liu 44 , K. Y. Liu 23 , Kai Liu 6 , P. L. Liu 29 , Q. Liu 6 , S. B. Liu 40 , 
X. Liu 22 , X. H. Liu 1 , Y. B. Liu 26 , Z. A. Liu 1 , Zhiqiang Liu 1 , Zhiqing Liu 1 , H. Loehner 21 , 
G. R. Lu 12 , H. J. Lu 14 , J. G. Lu 1 , Q. W. Lu 30 , X. R. Lu 6 , Y. P. Lu 1 , C. L. Luo 24 , M. X. Luo 45 , 
T. Luo 37 , X. L. Luo 1 , M. Lv\ C. L. Ma 6 , F. C. Ma 23 , H. L. Ma 1 , Q. M. Ma 1 , S. Ma 1 , T. Ma 1 , 
X. Y. Ma 1 , Y. Ma 11 , F. E. Maas 11 , M. Maggiora 43j4 ' 43C , Q. A. Malik 42 , Y. J. Mao 27 , Z. P. Mao 1 
J. G. Messchendorp 21 , J. Min 1 , T. J. Min 1 , R. E. Mitchell 17 , X. H. Mo 1 , C. Morales 
Morales 11 , C. Motzko 2 , N. Yu. Muchnoi 5 , H. Muramatsu 39 , Y. Nefedov 20 , C. Nicholson 6 , 
I. B. Nikolaev 5 , Z. Ning 1 , S. L. Olsen 28 , Q. Ouyang 1 , S. Pacetti 18B , J. W. Park 28 , M. Pelizaeus 37 
H. P. Peng 40 , K. Peters 7 , J. L. Ping 24 , R. G. Ping 1 , R. Poling 38 , E. Prencipe 19 , M. Qi 25 , 
S. Qian 1 , C. F. Qiao 6 , X. S. Qin 1 , Y. Qin 27 , Z. H. Qin 1 , J. F. Qiu 1 , K. H. Rashid 42 , G. Rong 1 , 
X. D. Ruan 9 , A. Sarantsev 20 > c , B. D. Schaefer 17 , J. Schulze 2 , M. Shao 40 , C. P. Shen 37 > d , 
X. Y. Shen 1 , H. Y. Sheng 1 , M. R. Shepherd 17 , X. Y. Song 1 , S. Spataro 43j4 ' 43C , B. Spruck 36 , 
D. H. Sun 1 , G. X. Sun 1 , J. F. Sun 12 , S. S. Sun 1 , Y. J. Sun 40 , Y. Z. Sun 1 , Z. J. Sun 1 , Z. T. Sun 40 
C. J. Tang 31 , X. Tang 1 , I. Tapan 35C , E. H. Thorndike 39 , D. Toth 38 , M. Ullrich 36 , G. S. Varner 37 

B. Wang 9 , B. Q. Wang 27 , K. Wang 1 , L. L. Wang 4 , L. S. Wang 1 , M. Wang 29 , P. Wang 1 , 

P. L. Wang 1 , Q. Wang 1 , Q. J. Wang 1 , S. G. Wang 27 , X. L. Wang 40 , Y. D. Wang 40 , Y. F. Wang 1 
Y. Q. Wang 29 , Z. Wang 1 , Z. G. Wang 1 , Z. Y. Wang 1 , D. H. Wei 8 , P. Weidenkaff 19 , 
Q. G. Wen 40 , S. P. Wen 1 , M. Werner 36 , U. Wiedner 2 , L. H. Wu 1 , N. Wu 1 , S. X. Wu 40 , 
W. Wu 26 , Z. Wu 1 , L. G. Xia 34 , Z. J. Xiao 24 , Y. G. Xie\ Q. L. Xiu 1 , G. F. Xu\ G. M. Xu 27 , 

H. Xu 1 , Q. J. Xu 10 , X. P. Xu 32 , Z. R. Xu 40 , F. Xue 15 , Z. Xue 1 , L. Yan 40 , W. B. Yan 40 , 
Y. H. Yan 16 , H. X. Yang 1 , Y. Yang 15 , Y. X. Yang 8 , H. Ye 1 , M. Ye 1 , M. H. Ye 4 , B. X. Yu 1 , 



1 



43 C. X. Yu 26 , J. S. Yu 22 , S. P. Yu 29 , C. Z. Yuan 1 , Y. Yuan 1 , A. A. Zafar 42 , A. Zallo 18A , Y. Zeng 16 , 

44 B. X. Zhang 1 , B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , H. H. Zhang 33 , H. Y. Zhang 1 , 

4 5 J. Q. Zhang 1 , J. W. Zhang 1 , J. Y. Zhang 1 , J. Z. Zhang 1 , S. H. Zhang 1 , X. J. Zhang 1 , 

46 X. Y. Zhang 29 , Y. Zhang 1 , Y. H. Zhang 1 , Y. S. Zhang 9 , Z. P. Zhang 40 , Z. Y. Zhang 44 , G. Zhao 1 , 

47 H. S. Zhao 1 , J. W. Zhao 1 , K. X. Zhao 24 , Lei Zhao 40 , Ling Zhao 1 , M. G. Zhao 26 , Q. Zhao 1 , 

48 S. J. Zhao 46 , T. C. Zhao 1 , X. H. Zhao 25 , Y. B. Zhao 1 , Z. G. Zhao 40 , A. Zhemchugov 20 ' a , 

49 B. Zheng 41 , J. P. Zheng 1 , Y. H. Zheng 6 , B. Zhong 1 , J. Zhong 2 , L. Zhou 1 , X. K. Zhou 6 , 

50 X. R. Zhou 40 , C. Zhu 1 , K. Zhu 1 , K. J. Zhu 1 , S. H. Zhu 1 , X. L. Zhu 34 , X. W. Zhu 1 , 

51 Y. C. Zhu 40 , Y. M. Zhu 26 , Y. S. Zhu 1 , Z. A. Zhu 1 , J. Zhuang 1 , B. S. Zou 1 , J. H. Zou 1 

52 (BESIII Collaboration) 

53 1 Institute of High Energy Physics, Beijing 100049, P. R. China 

54 2 Bochum Ruhr-University, Jf.Jj.180 Bochum, Germany 

55 3 Carnegie Mellon University, Pittsburgh, PA 15213, USA 

56 4 China Center of Advanced Science and Technology, Beijing 100190, P. R. China 

57 5 G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 

58 6 Graduate University of Chinese Academy of Sciences, Beijing 100049, P. R. China 

59 7 GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 

60 8 Guangxi Normal University, Guilin 541004, P. R. China 

61 9 GuangXi University, Nanning 530004, P.R.China 

62 10 Hangzhou Normal University, Hangzhou 310036, P. R. China 

63 11 Helmholtz Institute Mainz, J.J. Becherweg 4°~,D 55099 Mainz, Germany 

64 12 Henan Normal University, Xinxiang 453007, P. R. China 

65 13 Henan University of Science and Technology, Luoyang 471003, P. R. China 

66 14 Huangshan College, Huangshan 245000, P. R. China 

67 15 Huazhong Normal University, Wuhan 430079, P. R. China 

68 16 Hunan University, Changsha 410082, P. R. China 

69 17 Indiana University, Bloomington, Indiana 4^405, USA 

70 18 (A)INFN Laboratori Nazionali di Frascati, Frascati, Italy; 

71 (B)INFN and University of Perugia, 1-06100, Perugia, Italy 

72 19 Johannes Gutenberg University of Mainz, 

73 Johann-Joachim-Becher-Weg 45, 55099 Mainz, Germany 

74 20 Joint Institute for Nuclear Research, 141980 Dubna, Russia 

75 21 K VI/ University of Groningen, 9747 AA Groningen, The Netherlands 

76 22 Lanzhou University, Lanzhou 730000, P. R. China 

77 23 Liaoning University, Shenyang 110036, P. R. China 

78 24 Nanjing Normal University, Nanjing 210046, P. R. China 

79 25 Nanjing University, Nanjing 210093, P. R. China 
so 26 Nankai University, Tianjin 300071, P. R. China 
si 27 Peking University, Beijing 100871, P. R. China 

82 28 Seoul National University, Seoul, 151-747 Korea 

83 29 Shandong University, Jinan 250100, P. R. China 

84 30 Shanxi University, Taiyuan 030006, P. R. China 

85 31 Sichuan University, Chengdu 610064, P. R. China 

86 32 Soochow University, Suzhou 215006, P. R. China 

87 33 Sun Yat-Sen University, Guangzhou 510275, P. R. China 



2 



88 Tsinghua University, Beijing 100084, P- R- China 

89 3 5 (A) Ankara University, Ankara, Turkey; (B)Dogus University, 

90 Istanbul, Turkey; (C)Uludag University, Bursa, Turkey 

91 36 Universitaet Giessen, 35392 Giessen, Germany 

92 37 University of Hawaii, Honolulu, Hawaii 96822, USA 

93 38 University of Minnesota, Minneapolis, MN 55455, USA 

94 39 University of Rochester, Rochester, New York 14627, USA 

95 40 University of Science and Technology of China, Hefei 230026, P. R. China 

96 41 University of South China, Hengyang 421001, P. R. China 

97 42 University of the Punjab, Lahore- 54590, Pakistan 

98 43 (A)University of Turin, Turin, Italy; (B)University of 

99 Eastern Piedmont, Alessandria, Italy; (C)INFN, Turin, Italy 

100 44 Wuhan University, Wuhan 430072, P. R. China 

101 45 Zhejiang University, Hangzhou 310027, P. R. China 

102 46 Zhengzhou University, Zhengzhou 450001, P. R. China 

103 a also at the Moscow Institute of Physics and Technology, Moscow, Russia 

104 b on leave from the Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

105 c also at the PNPI, Gatchina, Russia 

106 d now at Nagoya University, Nagoya, Japan 

Abstract 



Hadronic transitions of Xcj n c ir + ir~ (J = 0, 1, 2) are searched for using a sample of 1.06 x 
10 s -0(3686) events collected with the BESIII detector at the BEPCII storage ring. The r] c is 
reconstructed with KgK^ir^ and K + K~tt° final states. No signals are observed in any of the three 
Xcj states in either n c decay mode. At the 90% confidence level, the upper limits are determined 
to be B(xc0 -»• f? c 7r + < 0.07%, B(xa -»• »fcT + < 0.32%, and B(x C 2 -> ?? c vr + 7r~) < 0.54%. 
The upper limit of B(xd ~^ f] c -K + Ti _ ) is lower than the existing theoretical prediction by almost an 
order of magnitude. The branching fractions of Xcj i^^i^ =t vr =F 7r + 7r _ , K + K~ 7r + 7r~7r° , luK + K~ 
and (pTT + TT~TT° (J = 0, 1, 2) are measured for the first time. 

107 PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Pq 
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108 I. INTRODUCTION 



109 Heavy quarkonia, both cc and bb bound states, have provided good laboratories for the 

no study of the strong interaction [l], 0] . For the hadronic transitions between the heavy quarko- 

iii nium states, Yan |3[ characterized it as the emission of two soft gluons from the heavy quarks 

n2 and the conversion of gluons into light hadrons. Based on this scheme, a series of decay rates, 

in such as the El-El hadronic transition, El-Ml hadronic transition, Ml-Ml hadronic transi- 

H4 tion have been calculated 0, @| • It has been shown that the multipole expansion can make 

us quite successful predictions for many hadronic transitions between the heavy quarkonia [H,0]. 

lie However, most of these studies are for the transitions among the 3 S± states; the hadronic 

in transitions of 3 Pj states are seldom explored. Using a sample of T(3S), CLEO measured 

us for the first time the transition rate of P-wave bottomonium Yf>j(2P) — >■ X5j(lP)/T7r 

n9 and the results are consistent with the theoretical predictions [5j. For the hadronic tran- 

120 sition of the P-wave charmonium states, there is only an upper limit of 2.2% at the 90% 

121 confidence level (C.L.) on the Xc2 — ► ti c tt + 7t^ transition rate recently reported by the BaBar 

122 experiment 0. The most promising process Xci ~^ Vc^^, which is dominated by an El- 

123 Ml transition, is calculated in the multipole expansion formalism, and a transition rate 

124 B(xd -> r) c irir) = (2.72 ± 0.39)% is predicted @. 

125 In this article, we search for Xcj — > r] c 7r + 7r~ with rj c decays into K s K ± ir T and K + K~tt°, 

126 where the K s is reconstructed in tt + ti^ and 7r° in 77 final states. We also report the first 

127 measurement of the branching fractions of Xcj — > K s K ± 'k t 7t + 7t~ , K + K~7T + 7T~7r°, uK + K~, 

128 and (j)n + ix~n°. 



129 II. THE EXPERIMENT AND DATA SETS 
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The data sample for this analysis consists of 1.06 x 10 8 events produced at the peak of 

1 of data were collected at a center-of-mass 



the ^(3686) resonance (10J; an additional 42 pb 
energy of ^=3.65 GeV to determine non-resonant continuum background contributions. 
The data are accumulated with the BESIII detector operated at the BEPCII e + e~ collider. 



The BESIII detector, described in detail in Ref. [Uj, has an effective geometrical accep- 
tance of 93% of Air. It contains a small cell helium-based main drift chamber (MDC) which 
provides momentum measurements of charged particles; a time-of-flight system (TOF) based 
on plastic scintillator which helps to identify charged particles; an electromagnetic calorime- 
ter (EMC) made of Csl (Tl) crystals which is used to measure the energies of photons 
and provide trigger signals; and a muon system (MUC) made of Resistive Plate Chambers 
(RPC). The resolution of the charged particles is 0.5% at 1 GeV/c in a 1 Tesla magnetic 
field. The energy loss (dE/dx) measurement provided by the MDC has a resolution better 
than 6% for electrons from Bhabha scattering. The photon energy resolution can reach 2.5% 
(5%) at 1 GeV in the barrel (endcaps) of the EMC. And the time resolution of TOF is 80 ps 
in the barrel and 110 ps in the endcaps. 

Monte Carlo (MC) simulated events are used to determine the detection efficiency, op- 
timize the selection criteria, and study the possible backgrounds. The simulation of the 
BESIII detector is geant4 12] based, where the interactions of the particles with the de- 
tector material are simulated. The ^(3686) resonance is produced with kkmc 13|, while 
the subsequent decays are generated with EvtGen The study of the background is 
based on a sample of 10 8 ^(3686) inclusive decays which are generated with known branch- 



ing fractions taken from the Particle Data Group (PDG) [15], or with lundcharm [16J for 
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152 the unmeasured decays. 



153 III. EVENT SELECTION 



154 A charged track should have good quality in the track fitting and be within the an- 

155 gle coverage of the MDC, |cos#| < 0.93. A good charged track (excludes those from K s 

156 decays) is required to be within 1 cm of the e + e~ annihilation interaction point (IP) trans- 

157 verse to the beam line and within 10 cm of the IP along the beam axis. Charged-particle 
i5B identification (PID) is based on combining the dE/dx and TOF information in the variable 

159 XpIdW = ( dE / dx —^-dE/dx mpcctad y + ^ TOF mcaaur ^-TOF cxpcctcd ^ 2 _ ^ ^2^) &nd the 

160 corresponding confidence levels Probpro^) are calculated for each charged track for each 
lei particle hypothesis i (pion, kaon, or proton). 

162 Photons are reconstructed from isolated showers in the EMC which are at least 20 degrees 

163 away from any of the charged tracks. In order to improve the reconstruction efficiency and 

164 the energy resolution, the energy deposited in the nearby TOF counter is included. Photon 

165 candidates are required to have the energy greater than 25 MeV in the EMC barrel region 

166 (| cos^l < 0.8), while in the EMC endcap region (0.86 < | cos#| < 0.92), the energy threshold 

167 requirement is increased to 50 MeV. EMC timing requirements are used to suppress noise 
lea and energy deposits unrelated to the event. 

169 K s candidates are reconstructed from secondary vertex fits to all the charged-track pairs 

170 in an event (assume the tracks to be tt). The combination with the best fit quality is kept, and 

171 the Kg candidate must have an invariant mass within 7 MeV/c 2 of the Kg nominal mass and 

172 the secondary vertex be at least 0.5 cm away from the IP. The reconstructed K s information 

173 is used as input for the subsequent kinematic fit. The ir° candidates are reconstructed from 

174 pairs of photons with an invariant mass in the range 0.120 GeV/c 2 < M( / -f-f) < 0.145 GeV/c 2 . 

175 In selecting ip(3Q86) —> / -fK s K ± 7r T 'R + 'ir~ , a candidate event should have at least six 

176 charged tracks and at least one good photon. After K s selection, the event should have 

177 exactly four additional good charged tracks with zero net charge. While in selecting 
i7B ^(3686) — > / -fK + K~7T + 7T~7r°, a candidate event should have four good charged tracks with 
179 zero net charge and at least three good photons. The r yK s K ± 7r T 7T + 7r~ [ r yK + K^^-k^-k ) 
iso candidate is then subjected to a four-constraint (4C) kinematic fit to reduce background 

181 and improve the mass resolution. Determination of the species of the final state parti- 

182 cles and selection of the best photons when additional photons (and 7r° candidates) are 

183 found in an event are achieved by selecting the combination with the minimum value of 

184 x 2 — xlc + Sj=i XpidC/)' where xlc * s the chi-square from the 4C kinematic fit. Events 

185 with xlc < 50 are kept as •yKgK ± TT T ir + n~ (•yK + K~7i + n~'K ) candidates. 

186 There is substantial background from ^(3686) —> X + J/ip decays. The ?/>(3686) — > 

187 7i + 7i^J/tp events are removed by requiring the recoil mass of any 7r + 7r~ pair to be outside 
las of a ±3<7 window around the J/ip nominal mass, where a is the resolution of the tt + tt~ 
189 recoil mass. The ^(3686) — > r/J/ip, r\ — >• 77r + 7r~ events are rejected if 0.535 GeV/c 2 < 
wo M(77r + 7r~) < 0.555 GeV/c 2 . To reject background with one more photon than the signal 

191 events {^^K s K ± y;^'K Jr 'K~ or , ~fyK + K~7T + 7T~7T ), x\q > 10 is required when a 4C kine- 

192 matic fit to ^KgK^ix^-n^-n^ or 77i^ + i^~7r + 7r _ 7r° is applied to the event; and to suppress 

193 ^(3686) ttVJ/V, J/i[) K + K-n + 7r~ background, M(K + K~'k + tx') > 3.125 GeV/c 2 

194 or M(K+K~7i+7i~) < 3.095 GeV/c 2 is required. 
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195 IV. DATA ANALYSIS 

we A. XcJ — > KgK ± 7r^ 1 7r + 7v~ and XcJ — > K+K _ 7r+7r _ 7T 

197 After the above selection, the invariant mass distributions of the hadron system are shown 

198 in Fig. [U and clear Xcj signals are observed with very low background level in the two modes. 




FIG. 1: Invariant mass spectrum of KgK^Tr^TT^ir - (left panel) and K + K~ "7r + tt~ it (right panel), 
together with the best fit results. The points with error bars are data, and the solid lines are the 
total fit results. The Xco, Xci, and Xc2 signals are shown as dotted lines, dash-dotted lines, and 
long-dashed lines, respectively; the backgrounds are in dashed lines. 



199 Using the inclusive MC events sample, the potential backgrounds from the 0(3686) 
decays which may contaminate Xcj — > KgK ± ii T 7i + 7i^ (x c j ~^ K + K~ n + n~ n ) are es- 
timated. Events from 0(3686) -> 7Xcj, Xcj -> KgKgTr + 7r~ and Xco -> KgK° s K + K~ 
(■0(3686) — > 'jXcj, Xcj KgK ± 7i T 7i°) may peak at the signal region, and the contributions 
203 from these peaking backgrounds are estimated using the detection efficiencies determined 
from the MC simulations and the corresponding branching fractions from previous measure- 
ments lH|, and then subtracted. These background events after final event selection are 



200 
201 
202 



204 
205 



206 listed in Table |T] and [TTJ the errors in the numbers of events are from the uncertainties of the 

207 detection efficiency and branching fractions. The other backgrounds are composed of dozens 

208 of decay modes and smoothly distributed in the full mass region (3.3 GeV/c 2 ~ 3.6 GeV/c 2 ); 

209 this kind of background is described by a second-order Chebyshev function. 



TABLE I: The number of remanent peaking background events (N^f 1 ) in Xc.i —> KgK^ir^ 



IT ' IT 



after final event selection. The branching fractions (B) are taken from PDG 15] . 



Decay modes 


ArPcak 
iV bk K 


sua 




Xco -> K%K%v + n- 


4.9 ±1.4 


(5.8 ± 1.1) x 10" 


-3 


Xd -> K° s K° s Tr+7r- 


0.1 ±0.1 


(7.2 ±3.1) x 10" 


-4 


Xc2 -> K° S K° S 7T+7T- 


0.6 ±0.3 


(2.4 ±0.6) x 10" 


-3 


Xco K° S K° S K+K- 


43.3 ± 15.8 


(1.4 ±0.5) x 10" 


-3 
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TABLE II: The number of remanent peaking background events (A^*) in XcJ_~^ K + K ir + ir tt' 



after final event selection. The branching fractions (£>) are taken from PDG 15] 



Decay modes 


ArPcak 
iV bk K 




[15J 


XcO -> K° S K^7T° 


123.8 ± 19.4 


±x( 


2.52 ± 0.34) x 10" 2 


Xcl -> K^K^n 


39.0 ±9.1 


Ix 


(9.0 ± 1.5) x 10~ 3 


Xc2 -> A° A^vr 


52.5 ±9.3 


Ix( 


1.51 ±0.22) x 10~ 2 



Data taken at y/s = 3.65 GeV are used to estimate backgrounds from the continuum 
process e + e~ — > qq. This kind of background is found to be small and uniformly distributed 
in the full mass region of interest in both decay modes, so the contribution can be represented 
by the smooth background term in the fit. 

In the measurement of B(x c j ~± K + K~ir + n~n G ), the contributions from intermediate 
states with narrow resonances such as 77, oj, and are excluded; the identification of such 
resonances is similar to that used in Refs. (l9l liol ] . The branching fractions of these decays 
Xcj — > r]K + K~ , Xcj — > uK + K~, Xcj — > 4>K + K~, and Xcj — > 4>Ti + ix~7t are measured using 
the same data sample. The fits to the invariant mass spectrum of 7r + 7r~7r° and K + K~ in 
the three Xcj signal regions (as defined in Sect. HVBj) ) are shown in Fig. [21 and the results 
are listed in Table IIHI The first errors are statistical and the second ones are systematic. 
The sources of the systematic errors are similar to those in the measurement of B(x c j 
K + K~7T + 7T~'n- ), as will be shown in Sect. [V] The branching fractions of Xcj ~^ ojK + K~ and 
Xcj — > 07r + 7T~7r are the first measurement, and those of other modes from this measurement 



are consistent within errors with the known PDG values [15] when available. Events from 
these decay modes are removed by requiring Im^+^-^o — m v \ > 15 MeV/c 2 , \m ir +^-^o—m u \ > 
40MeV/c 2 , \m 1T + 1T - 1T o—m ( j > \ > 15 MeV/c 2 and \m,K+K- — m <p\ > 15MeV/c 2 . The contribution 
from Xcj — > <P4>, <P ~^ K + K~, <p — > tt + 7t~7t° events is also removed by these requirements. 
The expected remaining events from these decay channels are also listed in Table IIHI and 
will be subtracted from the signal yields from the best fits. 

An unbinned maximum likelihood fit is applied to the invariant mass spectrum of 
KgK ± ir T 7i + 7i^ (K + K~7t + 7t~7t°) to extract the numbers of Xcj events in Fig. [TJ The Xcj 
signals are described by the corresponding MC simulated signal shape convolved with a 
Gaussian function G(fi, a) to take into account the difference in the mass scale and the 
mass resolution between data and MC simulation. The means (fi) and the standard devi- 
ations (cr) of the Gaussian functions are floated parameters in the fit. In the generation 
of the Xco MC events, the El radiative transition factor is included, where E y is the 
energy of the radiative photon in the ^(3686) rest frame. To damp the diverging tail due 
to the E^ dependence, a damping function E Eo+ ^ _ Eo )2 used by KEDR [l8| is introduced, 

where E = (fn^ 3686 ^ — m 2 c0 )/2m^(3 686 ). The backgrounds are described by a second-order 
Chebyshev function in both decay modes. 

The fit to the invariant mass spectrum of KgK ± TT T 7i + 7i^ yields 2837 ±64, 5180 ±75, and 
4560 ±71 signal events for Xco, Xa, and Xc2, respectively; while the fit to the K + K~7v + n^n 
modes yields 9372±130, 12415±126, and 11366±123 events for Xco, Xci, and Xc2, respectively. 
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FIG. 2: Fits to the invariant mass spectrum of 7r + 7r~7r (left panels) and K + K~ (right panels) 
in three Xcj mass region after Xcj 4"fti 4> ~^ K + K~, <ft — > 7r + 7r~7r° were rejected. From top to 
bottom are in XcO mass region, Xci mass region, and Xc2 mass region, respectively. 



244 B. XcJ — ► r? c 7r + 7r 

245 To study the Xcj rj c 7r + ir~ transitions, we define the signal regions for x c0 , Xci, and Xc2 as 

246 3.38 GeV/c 2 < M(KKinnr) < 3.45 GeV/c 2 , 3.48 GeV/c 2 < M{KKhtxtx) < 3.54 GeV/c 2 , 

247 and 3.54 GeV/c 2 < M{KKixixtx) < 3.60 GeV/c 2 , respectively. Here, KKimii is either 

248 KgK ± 7T T ii + ii^ or K + K~ n + it~ n Q . The invariant mass spectra of KgK ± n T (K + K~ir ) with 

249 KgK ± 7r T ii + ii^ (K + K~7T + 7T~7r°) in the three Xcj signal regions are shown in Fig. [31 In both 

250 decay modes, there are no significant r\ c signal in Xco and Xci decay; the r\ c signal observed 

251 in Xc2 decays is found to be mainly from ?/;(3686) — > ix + tx~ J J /ip — > 7?7 C , r] c — > KgK ± n T 

252 (r] c — > K + K~7i°), as described below. 

253 The potential backgrounds from ?/;(3686) decays are investigated with the inclusive 
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TABLE III: The number of background events (N^^ ) remained in Xcj K + K~7r + ir~7r° from 
modes with narrow intermediate states. The branching fraction of Xcj 4"^ ls taken from the 
BESIII measurement 19[] with both statistical and systematic errors; the other branching fractions 
are measured in this analysis. We also list PDG values 15[ in the last column for comparison. 



Decay modes 




Branching fraction 


r\JKj [ioj value 


XcO -> 


r]K + K~ 





forbidden by J p -conservation 


< 2.3 x 10 


Xcl -> 


T]K + K~ 


o.ol ± (J. 


/ o /loin z f 7\ w in — 4 
(3.48 ± 0.5 J) X 1U 


[6.6 ± l.U ) X 1U 


Xc2 -> 


r)K + K" 


1.9U ± 0.53 


(1.69 ± 0.45) x 10 


< 3.5 x 10 


XcO ->■ 


loK + K~ 


i of 1 o t i i n O/i 

13o.87 ± 10.84 


(1.94 ± 0.06 ± 0.20) x 10 




Xcl ->■ 


ujK+K- 


oy.z4 ± o.zu 


( ( .5/ ± U.oD ± U.04) X 1U 




Xc2 -> 


uK + K- 


57.26 ± 5.04 


(7.32 ±0.39 ±0.78) x 10" 4 




XcO ->■ 


(f>K+K- 


30.37 ±4.91 


(1.15 ±0.17) x 10~ 3 


(0.98 ±0.25) x 10~ 3 


Xcl ->■ 


(f>K+K- 


17.61 ±3.79 


(6.74 ±1.37) x 10~ 4 


(4.3 ± 1.6) x 10~ 4 


Xc2 -> 


(j)K + K- 


30.39 ± 5.27 


(1.14 ±0.16) x 10~ 3 


(1.55 ±0.33) x 10- 3 


Xco ->■ 




49.57 ±2.68 


(1.18 ±0.07 ±0.13) x 10~ 3 




Xcl ->■ 




28.59 ± 1.99 


(7.54 ±0.53 ±0.80) x 10~ 4 




Xc2 -> 


(f)lT + 1T~1T 


34.45 ± 2.36 


(9.25 ±0.63 ±0.97) x 10~ 4 




XcO 




0.97 ±0.17 


(8.0 ±0.3 ±0.8) x 10" 4 


(8.2 ±0.8) x 10" 4 


Xcl 




0.38 ± 0.08 


(4.4 ±0.3 ±0.5) x 10~ 4 




Xc2 




1.03 ±0.20 


(10.7 ± 0.3 ± 1.2) x 10" 4 


(11.4 ± 1.2) x 10~ 4 



254 MC events. The dominant backgrounds are the irreducible ^(3686) — > 'JXcj, Xcj 

255 K s K ± n T -K + it~ (xcj K + K~n + 'K~'K ) events. These events have the same final states 

256 as the signal events but KgK ± n T (K + K~ir°) are not from the decay of r\ c . In the Xc2 signal 

257 region, there are peaking backgrounds from the decay ^(3686) —> tt + tt~ J/if>, J /if; — > / ~fr] c , 

258 f] c — > KgK ± n T (r] c — > K + K~ir°). The energy of the transition photon in this decay is close 

259 to the energy of the transition photon in ^(3686) — > 7Xc2, an d the final states are the same 

260 as the signal events. The other backgrounds are composed of dozens of decay channels, each 

261 with a small contribution. The dominant backgrounds and the other backgrounds contribute 

262 a smooth component in the t] c mass region (2.7 GeV/c 2 ~ 3.2 GeV/c 2 ), so these backgrounds 

263 are described by a second-order Chebyshev function as shown in Fig. [3j In the Xc2 case, the 

264 peaking background has the same final state and similar kinematics as the signal events, so 

265 we use the same line-shape to describe both of them. Data taken at ^fs = 3.65 GeV are 

266 used to estimate backgrounds from the continuum process e + e~ — > qq. It is found that this 

267 background is small and uniformly distributed in the full mass region in both decay modes, 

268 so the contribution is neglected. 

269 An unbinned maximum likelihood fit is applied to the invariant mass spectrum of 

270 KgK ± 7T T (K + K~tt°) to extract the number of rj c events, as shown in Fig. [3j The r] c signal is 

271 described by a MC simulated line shape with the detector resolution included, and the reso- 



272 nance parameters of i] c are fixed to the latest measurement from the BESIII experiment [17 

273 The background (except the peaking background in the Xc2 signal region) is described with 

274 a second-order Chebyshev polynomial function in both decay modes in the three Xcj signal 

275 regions. 
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FIG. 3: Invariant mass spectra of KgK^-ir^ (top row) and K + K~ir° (bottom row) with KK-k-rtt 
in XcO (left panel), Xci (middle panel), and Xc2 (right panel) signal regions and the fit results. Dots 
with error bars are data; the solid lines are the total results from the best fits to the invariant mass 
spectrum. The r/ c signals are shown in dash-dotted lines (in Xc2 mass region, the contribution from 
the peaking background is not removed.); the backgrounds as dashed lines. 



276 As there is no significant r\ c signal in any of the three Xcj states in either rj c decay mode, 

277 we set upper limits on B(x c j ?7c 7r+7r_ ) using the probability density function (PDF) for 
27B the expected number of signal events. In the Xco and Xci signal regions, the likelihood 

279 distributions in the fitting of the invariant mass spectra in Fig. [3] are taken as the PDFs 

280 directly. They are obtained by setting the number of r] c signal events from zero up to a 

281 very large number. In the Xc2 signal region, the likelihood distribution also contains the 

282 contribution from the peaking background. Using the known branching fractions [l5|, the 

283 detection efficiency from MC simulation, and the number of ^(3686) events, the expected 

284 peaking background are 45. 7± 11.6 in KgK ± ir T and 34.4±8.7 in K + K~ir°. Here, the errors 

285 include the uncertainties in the detection efficiency and the branching fractions. Then the 

286 PDF of signal is extracted with the PDF of the peaking background (Gaussian distribution 

287 with mean set to the expected number of peaking backgrounds, sigma set to its error) and 

288 the PDF from the fit. The systematic uncertainties are considered by smearing the PDF in 

289 each decay with a Gaussian. The upper limit on the number of events at the 90% C.L. is 

290 defined as N up , corresponding to the number of events at 90% of the integral of the smeared 

291 PDF. In each decay mode in the three Xcj states, the fit-related systematic errors on the 

292 number of signal yield are estimated by using different fit ranges, different orders of the 

293 background polynomial, and different r] c line shapes with the parameters of r] c changed by 

294 one standard deviation [l7j]; the maximum N np is used in the upper limit calculation. 
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29 5 V. SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties in the measurement of B(x c j — > KKtttttt) and B(x c j — > 
r] c ir + TT~) are summarized in Tables HVl and N] respectively. The systematic errors related to 
the MDC tracking (2% per track for those from IP), photon reconstruction (1% per photon), 



and 7T° reconstruction (1%) are estimated with control samples [20|, [21[; the errors in the 
branching fractions of ^(3686) — > ^Xcj and t] c — y KKir are taken from the PDG [l5[ and 
are propagated to the Xcj branching fraction measurement; and a 2% uncertainty is taken 
for each decay due to the limited statistics of the MC samples used. There is an overall 4% 
uncertainty in the branching fraction associated with the determination of the number of 
■0(3686) events in our data sample (Io| . 



TABLE IV: Systematic errors (in %) in B(x c j -> K^K^tt+it-) and B(x c j -> K + K~ -k + -k~ tt q ) . 



Sources 


KsKVir+w- 

XcO Xcl Xc2 


K + K~ir+ 

XcO Xcl 


7T 7T° 

Xc2 


MDC tracking 




8.0 








8.0 






Photon reconstruction 




1.0 








3.0 






MC statistics 


1.1 


1.4 


1 


5 


1.0 


1.4 


1 


4 


Kg reconstruction 


1.4 


1.6 


1 


7 










7T° reconstruction 












1.0 






Kinematic fit 


1.5 


1.9 


1 


7 


0.4 


0.4 





2 


Damping function 


0.5 


0.1 





1 


0.4 


0.1 





1 


Intermediate states 


1.0 


1.0 


1 





4.0 


4.0 


4 





Fitting range 


1.0 


0.4 





2 


0.4 


0.4 





7 


Background shape 


1.4 


0.7 





6 


1.3 


0.7 





4 


£(0(3686) -> 1XcJ ) 


3.2 


4.4 


3 


9 


3.2 


4.4 


3 


9 


Number of '0(3686) events 




4.0 








4.0 






Total 


10.1 


10.5 


10.3 


10.9 


11.3 


11.2 



305 A. Kg reconstruction 

306 The uncertainty in the Kg reconstruction arises from three parts: the geometric accep- 

307 tance, the tracking efficiency, and the efficiency of K s selection. The first part is estimated 

308 using MC simulation, and the other two are studied using J/ip — y K* ± K T , K* ± — y K ^. 

309 By selecting a pair of i^ ± 7r =F , the recoil mass spectrum shows a clear K° signal. The effi- 

310 ciency of K s reconstruction is calculated with B ( K o ^, n + n -] x ( ni+n2 y2 > wri ere n\ is the number 
3n of K° obtained from a fit to the K ± ir T recoiling mass when there is a K s reconstructed in 

312 the recoil side that satisfies the K s selection, and is the number of K° from fitting to the 

313 K ± ir T recoiling mass spectrum when no K s candidate satisfies the K s selection. The differ- 

314 ence in the efficiency of K s reconstruction ( £ data / £ MC — 1) as a function of K s momentum is 

315 shown in Fig. HI The difference in the K s reconstruction between data and MC simulation 

316 is fitted with a linear function of the K s momentum as shown in Fig. H] together with the 
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TABLE V: Systematic errors (in %) in B{x c j Tj c ir + 'K ) in r/ c — > K S K -k^ and rj c — > K + K n { 
decay modes. 



Sources 


?7c ~~ ' 

XcO 


Xcl Xc2 


Vc 

\ 


■ K + K-ir° 

Xcl Xc2 






8.0 






8.0 




P n ntnn rprnn t^yi i inn 

J. HVJuVJll 1 V-^UllQ ul U^ulUll 




1.0 






3.0 




MC statistics 


1.8 


1.5 


1.6 


1.8 


1.5 


1.7 


K s reconstruction 


2.4 


2.3 


2.2 








7T° reconstruction 










1.0 




Kinematic fit 


1.5 


1.7 


1.8 


0.9 


0.2 


0.4 


£(^(3686) 7Xcj) 


3.2 


4.4 


3.9 


3.2 


4.4 


3.9 


£(r/ c -»• IOTtt) 




8.4 






8.4 




Number of -0(3686) events 




4.0 






4.0 




Total 


13.2 13.5 


13.4 


13.2 13.5 


13.4 



317 ±lcr envelops. Since the difference between data and MC is significant, we do a correction 

318 to the signal MC according to the momentum of K s , and the uncertainty of this correction 

319 is taken as the systematic error. 

320 The systematic errors in B(xcj — > K s K ± ix t -k + ii~) are found to be 1.4%, 1.6%, and 1.7%, 

321 for J = 0, 1, and 2, respectively; while for B(xcj 7 7c 7r+7r ~) ) they are 2.4%, 2.3%, and 2.2% 

322 for Xco, Xci, and Xc2, respectively. 




P(K°)(GeV) 

FIG. 4: The difference in the K s reconstruction efficiency between data and MC simulation (points 
with error bars), together with the fit to the difference with a linear function of momentum. The 
solid line is from the best fit and the dashed lines are the ±la envelopes of the best fit. 



323 B. Kinematic fit 

324 In the MC simulation, the model is much simpler than the real detector performance, and 

325 this results in differences between data and MC simulation in the track parameters of photons 

326 and charged tracks. The simulation of the photon has been checked in another analysis (Io| . 
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327 which shows good agreement between data and MC simulation. For the charged tracks, 

328 careful comparisons with purely selected data samples indicate that the MC simulates the 

329 momentum and angular resolutions significantly better than those in data, while the error 

330 matrix elements agree well between data and MC simulation. This results in a much narrower 

331 x\c distribution in MC than in data, and introduces a bias in the efficiency estimation. We 

332 correct the track helix parameters of MC simulation to reduce the difference between data 

333 and MC simulation. 

334 We use J /if; — > 0/o(98O), 4> — > K + K~, / (980) —> 7i + n~ as a control sample to study 

335 the difference on the helix parameters of charged tracks between data and MC, as this 

336 channel has a large production rate, very low background, and has both pions and kaons. 

337 We find that the pull distributions of data are wider than MC simulation and the peak 

338 positions are shifted. These obvious differences between data and MC suggest wrong track 

339 parameters have been set in MC simulation. The helix parameters of each track in the MC 

340 simulation are enlarged by smearing with a Gaussian function G(p, a) = G((pf ata — /if 10 ) x 

341 Vu, ^/(af^/af 10 ) 2 — 1 x Vu), where % = {dp, O > K , dz,tg\} is the i-th helix parameter of 

342 the track and V is the corresponding covariance matrix. Here dp is the distance from the 

343 pivot to the orbit in the x-y plane, <fro is the azimuthal angle specifies the pivot with respect 

344 to the helix center, k is the reciprocal of the transverse momentum, dz is the distance of 

345 the helix from the pivot to the orbit in the z direction, and tgX is the slope of the track. 

346 The correction factors pf atSj , pf 10 , <rf ata , and af 10 are the means and resolutions of the pull 

347 distributions of the data and MC obtained from control samples. 

348 The correction factors are listed in Table IVII If the correction is perfect, xtc distribu- 

349 tions of data and MC simulation will be consistent with each other; however, from the 

350 comparisons of many final states, we find that the agreement between data and MC sim- 

351 ulation does improve significantly but differences still exist. This indicates that the effect 

352 is from multiple sources, and our procedure cannot solve all the problems. In our analy- 

353 sis, we take the efficiency from the track-parameter-corrected MC samples as the nominal 

354 value, and take half of the difference between MC samples before and after the correction 

355 as the systematic error from the kinematic fitting. This is a very conservative estimation. 

356 The comparison of xtc distributions between data and MC simulation before and after 

357 the track-parameter-correction are shown in Figs. and for Xcj ~^ KgK 7r^7r + 7r~ and 

358 K + K~-k + 7t~7t°, respectively. 

TABLE VI: Correction factors extracted from pull distributions using a control sample of J /if) — > 
<£/o(980), <f) -> K+K-, / (980) -> tt+tt'. 





00 


K 


tgX 




^data _ ^MC ^data /^MC 


^data _ ^MC ^.data y^-MC 


^data _ ^MC ^data y^MC 


K+ 


-0.04 ±0.03 1.19 ±0.02 


-0.24 ±0.03 1.28 ±0.02 


-0.38 ±0.01 1.25 ±0.02 


K~ 


0.06 ±0.03 1.21 ±0.02 


0.25 ±0.03 1.25 ±0.02 


-0.36 ±0.01 1.21 ±0.02 


7T+ 


-0.06 ±0.03 1.25 ±0.02 


-0.10 ±0.03 1.31 ±0.02 


-0.36 ±0.01 1.25 ±0.02 


7T~ 


-0.02 ±0.03 1.23 ±0.02 


0.10 ±0.03 1.27 ±0.02 


-0.36 ±0.01 1.21 ±0.02 



359 The systematic errors in B(xcj — > K^K ± n T it + ii ) (B(xcj K + K n + n 7T )) are 1.5%, 

360 1.9%, and 1.7% (0.4%, 0.4%, and 0.2%) for J = 0, 1, and 2, respectively; and 1.5%, 1.7%, 
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FIG. 5: Comparison of xlc between signal MC and data for -0(3686) — > 7Xcj, Xcj ~^ 
KgK^ir^^ir^ . The points with error bars are data, and the solid lines are MC simulation. 
Left panel: signal MC without track-parameter-correction; Right panel: signal MC after track- 
parameter-correction. From top to bottom are XcO, Xd, and Xc2, respectively. 



361 and 1.8% (0.9%, 0.2%, and 0.4%) systematic uncertainties are assigned to B(xcj t] c 7t + 'k ) 

362 with 7] c ->■ KgK ± 7i T (K + K~tt°) for J = 0, 1, and 2, respectively. 



363 C. Uncertainty from damping factor 

364 In the fit to the invariant mass spectrum of KgK ± TT T ir + n^ and K + K^7Y + n^7r , the 

365 damping function used by KEDR is adopted. Another damping factor used by CLEO [22| 

366 is e~ E ~<l^ with = (0.0650 ± 0.0025) GeV determined from their fit. Using this damping 

367 function with f3 = (0.097±0.024) GeV which is extracted from fitting Xco data, the differences 
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FIG. 6: Comparison of xlc between signal MC and data for -0(3686) — > 7Xcj> Xcj ~^ 
K + K~7t + 7T~7r° . The points with error bars are data, and the solid lines are signal MC plus possible 
background estimated using inclusive MC sample. Left panel: signal MC without track-parameter- 
correction; Right panel: signal MC after track-parameter-correction. From top to bottom are XcO> 
Xd, and Xc2, respectively. 



368 on the branching fractions of Xcj K s K ± ix t 'k + 'k~ (xcj K + K~7r + -K~7r°) are assigned to 

369 the systematic error due to damping function, which are 0.5%, 0.1%, and 0.1% (0.4%, 0.1%, 

370 and 0.1%) for J = 0, 1, and 2, respectively. The effect for Xci and Xc2 is small since the two 

371 states are very narrow. 
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372 D. Uncertainty from intermediate states 

373 The detection efficiencies for the measurement of B(x c j — > KgK ± TT T 7T + 7T~) and B(xcj — > 

374 K + K~ir + ff~7r°) are estimated using the MC simulation with Xcj decay to KgK ±/ K T 7T + 7T~ 

375 and K + K~7r + 7T~7T° generated according to pure phase space distribution. From the data, 

376 we see broad intermediate states such as K* and p in the invariant mass spectra of Kit 

377 and mi. The branching fractions of Xcj ~* KgK ± 7r T ir + ir~ (x c j K + K~ n + it ) via 

378 these intermediate states are measured by fitting the invariant mass spectra of Ktt and 

379 Tux. An alternative signal MC sample is generated with all possible intermediate states 

380 and corresponding branching fractions to determine the efficiency. The efficiency difference 

381 between this sample and the phase space sample is about 1.0% for Xcj KgK ± 7r T 7T + ir~ and 

382 4.0% for Xcj — > K + K~7r + 7T~7r°; these are taken as the systematic error due to intermediate 

383 states. 

384 E. Uncertainty from fitting 

385 The systematic uncertainty due to the fitting range is estimated by fitting the invariant 

386 mass spectrum in the range 3.25 GeV/c 2 ~ 3.61 GeV/c 2 and 3.35 GeV/c 2 ~ 3.60 GeV/c 2 . 

387 The biggest differences in the branching fractions are assigned as errors, which are 1.0%, 

388 0.4%, and 0.2% for Xco, Xci, and Xc2, respectively, in KgK ± 7i T 7Y + 7T^ decay; and 0.4%, 0.4%, 

389 and 0.7% for Xco, Xci, and Xc2, respectively, in K + K~ix + n~ix° decay. The background shape 

390 is changed from a second-order Chebyshev polynomial function to a third-order Chebyshev 

391 polynomial function, and the differences are taken to be the systematic errors, which are 

392 1.4%, 0.7%, and 0.6% for Xco, Xa, and Xc2, respectively, in KgK ± ir T 7r + n^; and 1.3%, 0.7%, 

393 and 0.4% for Xco, Xa, and Xc2 in K + K~7r + 7i~ir°. 

394 VI. RESULTS AND DISCUSSION 

395 Using the numbers of signal Xcj events from the fits, together with the corresponding 

396 efficiencies, the branching fractions of Xcj — > KgK ± 7i T 7i + 7i~ (x c j — > K + K^tx + -k~tx q ) are 

397 determined and listed in Table I VIII In the branching fractions of Xcj K + K~7T + 7T~7r°, 

398 contributions from narrow resonances rj, u, and <fi are subtracted. All these are first mea- 

399 surements, and the branching fractions are at the 1% level. Comparing the two decay 

400 modes, we found the ratio of the branching fractions is around one-half which may be a con- 

401 sequence of isospin symmetry. We also measured the branching fractions of Xcj ooK + K~ 

402 and Xcj 07T + 7r~7r° for the first time, the results are listed in Table [VTT] also. 

403 With the upper limit on the numbers of events at the 90% C.L. in Xcj f] c ^ + ^ ~, 

404 r] c — > KgK ± ir T (r] c — > K + K~7T°), as well as the corresponding efficiencies, the upper limits 

405 on the branching fraction of Xcj rj c 'K + 'iT~ in the two decay modes are determined, as 

406 listed in Table IVIIII We give a more stringent constraint on the Xc2 rj c ii + n~ branching 

407 fraction than BaBar does [8[. The theoretical prediction of B(xd Vc^ +7r ~) is also listed in 

408 Table rVTH[ which is larger than our measurements. We note that the theoretical prediction 

409 uses experimental results as input to normalize the parameters in the model. For example, 

410 the parameter a M loi E is extracted by comparing the branching fraction of ^(3686) — > h c 7i° 
4n between the theoretical calculation and the experimental measurement. This makes the 
412 prediction highly dependent on the former experimental results and theoretical models. 
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TABLE VII: The results for B(x c J -> K^K^ir+ir-), B(x c J ->• ir + K"7r + 7r-vr°), £(x c J -> 
wif+A' - ), and B(x c j — > (/>vr + 7r _ 7r°). The first errors are statistical and the second ones are sys- 
tematic. 



Decay mode 


^signal £ ( % ) £(xl0~ 3 ) 


XcO -> K^ttVtt- 
XcO ->• i4T + K~7r + 7r-7r 
XcO - ojK + K~ 

XcO -> 07T + 7r"7r° 


2789 ±66 9.30 4.22 ± 0.10 ± 0.43 
9031 ± 132 10.34 8.61 ± 0.13 ± 0.94 
1414 ±42 8.04 1.94 ±0.06 ±0.20 
538 ±29 9.16 1.18 ± 0.07 ± 0.13 


Xci -> A + A~7r + 7r~7r° 
Xci ->■ ujK + K~ 
Xci ->■ </>7r + 7r~7r° 


5180 ±75 10.21 7.52 ±0.11 ±0.79 
12256 ± 127 11.10 11.46 ± 0.12 ± 1.29 
628 ±29 9.34 0.78 ± 0.04 ± 0.08 
373 ±26 10.50 0.75 ± 0.06 ± 0.08 


Xc2 -»• #"^±7^+^- 
Xc2 - > A + A~7r + 7r~7r° 
Xc2 -> ujK + K~ 

Xc2 -> <M + 7r~7T° 


4559 ±71 9.76 7.30 ± 0.11 ± 0.75 
11189 ± 124 10.48 11.69 ± 0.13 ± 1.31 
512 ±27 8.58 0.73 ± 0.04 ± 0.08 
408 ±28 9.88 0.93 ± 0.06 ± 0.10 



TABLE VIII: Upper limits at the 90% C.L. on B{x c j — > r) c Tr + ir~) in the two r\ c decay modes. N &t 
is the number of events from the fits shown in Figj3j In Xc2 case, A fit includes the contribution 
from the peaking background ^(3686) — > ir + ir~ J/ip, J/tp — >■ 77/c r\ c — > KgK ± TT^ 1 (K + K~7r ). 



Decay mode 






e (%) B^(xcj 


-> r]cTr + ir-) (%) 


6 theory (XcJ ->• f?c^ + (%) 


XcO ~ 




0.0 ±4.6 


6.8 


6.29 


0.07 




XcO - 


-> (K + K-tt°)tt + tt- 


0± 15 


33.6 


6.82 


0.41 




Xci - 


-> (K^KV^tt+tt- 


18 ± 17 


48.7 


9.45 


0.32 


1.81 ±0.26 


Xci - 


-> Ik + K-tt )tt + tt~ 


6 ±25 


50.0 


9.82 


0.44 


Xc2 ~ 




77 ± 19 


64.1 


7.72 


0.54 




Xc2 " 


-> (K+K-Tr )7r + TT- 


89 ±26 


105.4 


7.83 


1.23 
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